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^ ■ Abstract 

Q \ I summarize recent results in lattice supersymmetry with special attention to 
^ ; = 1 Super Yang-Mills (SYM) theory. 
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Supersymmetric gauge theories exhibit fascinating phenomena which have attracted much attention 
since some results UJ about the ground state structure of some of these models have appeared. For 
, this reason, much more effort has been dedicated to formulating a lattice version of supersymmet- 

O \ ric theories (for some recent results, see for example, in the two dimensional Wess-Zumino model, 

in HI 4 1, four dimensional Wess-Zumino model |5|, four dimensional = 1 SYM theory using Wil- 
ff^ . son fermions [6.^,8,^ JO 1, chiral fermions [11 J and other supersymmetric theories, as for example, in 

O; C2l[l3l[l4l[T5l[l6l[l7l[i8]|. Unfortunately, I do not have time to cover all these subject and I refer the 

' reader to IT9ll . for a complete list of references. 

The dynamics of = 1 SYM is similar to QCD: confinement of the colored degrees of freedom 
and spontaneous chiral symmetry breaking. The mechanism of spontaneous chiral symmetry breaking is 
(~] , quite peculiar in this model: U{\)\i?< broken by the anomaly to Z2Nc (^c is the number of colors) and it 

is spontaneously broken down to Z2 due to a non-zero value of the gluino condensate l20ll . Concerning 
the low-lying mass spectrum, there are various scenarios in the literature (with different predictions) 
which can be directly investigated on the lattice in order to check their validity. The lattice is a fantastic 
\ laboratory which can help us to better understand whether supersymmetry is or not a symmetry of the 

Nature. Also the investigation of the Ward-Takahashi identity (WTi) is a non-perturbative effect and has 
been studied on the lattice both, numerically and in lattice perturbation theory. For the other way around, 
it is also interesting to study = 1 SYM far away from the supersymmetric limit in order to check the 
predictions of QCD for the glueball states, etc. Therefore, placing the QCD in a wider variety of theories 
as SYM can give us more information about it. 

The outline of my talk is the following. First, an overview of the physics of the continuum A^ = 1 
SYM theory together with some of the non-perturbative phenomena are presented. Then, the difficulties 
in formulating a lattice version of supersymmetric gauge theories are described. Results and comparison 
when using either Wilson fermions or domain wall fermions are given. 

2. A^ = 1 SYM CONTINUUM THEORY: NON-PERTURBATIVE RESULTS 

The continuum action density for A^ = 1 SYM with a gauge group SU{Nc) reads 

C = --F-(x) F-lx) + \ A»(x)7^(P^A(x))» + m^^ (x)A»(x) , (1) 
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where A" is a 4-component Majorana spinor that satisfies the Majorana condition A" = A"^C. The 
gluon fields are represented by = -igA^^T" and F^^ = -igF^^T". V^X" = d^X" + y/abc^jlA^ is 
the covariant derivative in the adjoint representation. 

For nig = this action has a supersymmetry respect to the continuum supersymmetric transfor- 
mations, 

5A^{x) = -2gX{x)-ff,e , 

1 

SXix) = — aprFpr{x)e, 
9 

— % 

5X{x) = -eaprFpr{x) , (2) 

g 

where apr = |[7p, 7t], A = A"T" and e is a global Grassmann parameter with Majorana properties. 

These transformations relate fermions with bosons, they leave the action invariant, which means 
that the variation of the lagrangian under these tranformations is proportional to a total derivative of a 
current (then, the variation of the action is zero). Also, they commute with the gauge transformations so 
that the resulting Noether current, 5^, is gauge invariant. For A'^ = 1 SYM theory the supercurrent is 
given by, 

Sp{x) = -F^,{x)aprj^X''{x) . (3) 

It is also expected that the supercurrent satisfy a renormalized supersymmetric WTi, dpSj} = 
2mjiXR, where XR = ^x^' ^ = ^-^pcrfpcrA" and mji is the renormalized gluino mass. We have 
supersymmetry for niR = 0, while the non-vanishing of mji describes a soft breaking of supersymmetry 
(this can be easily seen applying the supersymmetric transformations over the mass term in the density 
action: We see that it is not possible to write down this term as a total derivative of an operator). It is 
generally assumed that supersymmetry is not anomalous and that only the mass term is responsible for a 
soft breaking. In 1211 however, the question of whether non-perturbative effects may raise a supersym- 
metry anomaly has been investigated. Only a study of the continuum limit of the lattice supersymmetric 
WTi can shed some light on this question. I will be back to this point later on. 



2.1 Chiral symmetry breaking 

A non-zero gluino mass term in Eq. Q, Cmass = rrigX'^X"', breaks supersymmetry softly (that means 
that all the non-renormalization theorem and cancellation of divergencies are preserved). In the massless 
case, the global chiral symmetry is U{l)x, which is broken by the anomaly, 

o t5 _ pupa pa pa /^x 

~ 327r2 ^i-iy^P^ '^^^ 

(where = A7^75A), leaves a Z2Nc subgroup unbroken. This Z2Nc symmetry group is itself broken 
down to Z2 due to the non-zero value of the gluino condensate. The conseguence of this spontaneous 
chiral symmetry breaking is a first order phase transition at mg = 0. That means the existence of Nc 
degenerate ground states with different orientations of the gluino condensate (A; = 0,---,A^c — 1), 

\ o 27rifc 

A^A"" ) = cK^e— . (5) 



Here A is the dynamical scale of the theory which can be calculated on the lattice, for example, while c 
is a numerical constant which depends on the renormalization scheme used to compute A. Eq. ^ shows 
up the dependence on the gauge group. For SU (2), two degenerate ground states with opposite sign of 
the gluino condensate appear |^, while for 5*^7(3), we have three degenerate ground states at k = k, 
(for a first numerical study see 122^). 
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2.2 The value of the gluino condensate 

The value of the gluino condensate in the pure = 1 SYM theory has been calculated in the eighties 
by using two different methods. One is based on strong coupling instanton calculations. 123 1 , while the 
second one is based on weak coupling instanton calculations, [24s|. They give different result for c (in 
Eq. ^) and this was known as the | puzzle. Various discussion about the validity of both methods 
can be found in the literature [25 1. More recently, a third elegant method |26| calculates the gluino 
condensate directly in the semiclassical approximation. This method gives results in agreement with 
the weak coupling instanton approximation fTM [T/l and confirm the correctness of the weak coupling 
instanton result. 

2.3 Light hadron spectrum 

= 1 SYM theory describes interactions between gluons and gluinos. In analogy with QCD one expect 
that the spectrum of the model consists of colorless bound states of those fundamental excitations, namely 
glueballs (gg), gluinoballs (AA) and gluino-glueball (gX). These fields can combined to form the chiral 
supermultiplet, 

S{y) = 4>{y) + V26x{y) + e^F{y) , (6) 

where (p represent the scalar and pseudoscalar gluinoballs while x is their fermionic partner. F is an 
auxiliary field. It is tempting to think that F represent the scalar and pseudoscalar glueballs when treating 
it as a dynamical field in the minimal Veneziano and Yankielowicz (VY) action, but this cause paradoxes 
l28l . For = 1 SYM theory a low energy effective action has been proposed by VY l29l . 

£e// = -{S^S)'/% + 7[(51og4 -S)\f + h.c] . (7) 

Expanding the effective VY action around its mininum, it is found that the low energy spectrum forms a 
supermultiplet, consisting of a scalar meson A'^A" (the adjoint— /o), a pseudoscalar meson A'^75A" (the 
adjoint— 77/), and a spin-1/2 gluino-glueball particle, the x ■ Glueballs are absent in this formulation. 
In the supersymmetric limit these masses are degenerate, while the introduction of a rrig 7^ breaks 
supersymmetry softly and leads to a splitting of the multiplet. Generalization of the VY effective action 
in order to include glueballs are discussed, for example, in fSOl lBTIl . These generalizations give results 
which are in substantial agreement with lattice simulation using Wilson fermions li7ii32J. 

3. LATTICE FORMULATION OF = 1 SYM THEORY 

There are several important motivations for defining a field theory on the lattice. The most important 
one is to regularize a "perturbative divergent theory" in a non-perturbative way. However, formulating a 
supersymmetric field theory on the lattice seems to be a difficult task. The first difficulty has to do with 
the failure of the Leibniz rule |33|. The lattice version of the Leibniz rule is order 0(a), thus, breaks 
supersymmetry. In fact, if we discretize in a trivial way the Eqs. (1 1121) we see that the variation of the 
lagrangian under the discretized lattice transformations is not proportional to a total lattice derivative. 
Only in the continuum limit one can recover the continuum results plus terms order 0{a) (depending 
on the discretization prescription). This is not a severe problem due to the fact that terms that break 
supersymmetry are irrelevant, i.e., they cancel in the continuum limit and therefore there is no need for 
fine tuning in order to eliminate their contributions. 

Another problem are the scalar mass terms in the supersymmetric theory that breaks supersym- 
metry. Since these operators are relevant, fine tuning is necessary in order to cancel their contributions. 
In the case of A^ = 1 SYM theory only one scalar mass term is present which has to be fine tuning. But 
in the case of extended supersymmetries, there will be a plethora of scalar fields and the situation would 
become drammatic if one would have to do a fine tuning for each of these operators. 



The last problem is also well known in QCD. A naive regularization of fermions originates the 
doubling problem |34|. This is also problematic for supersymmetry not only because the discretization 
results in a wrong number of fermions (16 instead of 1), but also results in a violation of the balance 
between bosons and fermions as supersymmetry requires. This problem can be treated as in QCD. In the 
following, I summarize some nice numerical results using the Wilson formulation and the domain wall 
approach (for a detailed review on these subjects, see tT9ll9ll35ll 

3.1 N = 1 SYM theory with Wilson fermions 

In the Wilson formulation of Curci and Veneziano l36l . it is proposed to construct a non supersymmetric 
discretized = 1 SYM theory with a supersymmetric continuum limit. In this formulation, supersym- 
metry is broken by the lattice itself, by the Wilson term and a soft breaking due to the gluino mass is 
present. Supersymmetry is recovered in the contimuum limit by tuning the bare parameters g and the 
gluino mass rrig (through the hopping parameter) to the supersymmetric limit. The supersymmetric (and 
chiral) limit are both recovered simultaneously at rrig = 0. 

The strategy using Wilson fermions is the following. By studying the pattern of chiral symmetry 
breaking, through the study of the first order phase transition of the gluino condensate it is then possible 
to determine the value of the critical hopping parameter which correspond to the supersymmetric limit 
(rrig = 0). In |6 , 22 1, for a fixed value of (3, it is introduced a gluino mass term that breaks supersymmetry 
and then it is tuned in order recovered supersymmetry in the continuum limit. At the supersymmetric 
(chiral) limit, a first order phase transition (or a crossover) should show up as a two double peak structure 
in the distribution of some order parameter (the gluino condensate, in this case), indicating that the 
corresponding k(. is the critical hopping parameter corresponding to the supersymmetric limit. How do 
we know we are restoring supersymmetry in the continuum limit? 

• This can be achive for example, by investigating the low-lying mass spectrum and comparing with 
theoretical predictions. In IX 32 1, the unquenched simulations using the multi-bosonic algorithm 
proposed by Liischer l37l with a two-step variant (the TSMB) algorithm l38ll . near the value of kc, 
have been performed. An accurate study of this issue is non-trivial not only from the computational 
point of view but also due to several different theoretical formulations. An independent method 
for checking the supersymmetry restoration would be demanding. 

• Another independent way to study the supersymmetry restoration is through the study of the su- 
persymmetric WTi. This has been achieved both numerically (U and in lattice perturbation theory 

cniioi. 

3.2 N = 1 SYM theory with Domain Wall fermions 

= 1 SYM theory has been also studied on the lattice using the domain wall fermion (DWF) approach. 
Application of DWF in supersymmetric theories has been explored in l4ni . The DWF approach is defined 
by extending the space-time to five dimensions. Also a non-zero five dimensional mass or domain wall 
height mo, which controls the number of flavors, is present. The size of the fifth dimension, Lg, and 
free boundary conditions for the fermions are implemented. As a result, the two chiral components of 
the Dirac fermion are separated with one chirality bound exponentially on one wall and the other on the 
opposite wall. For any value of a the two chiralities mix only by an amount that decreases exponentially 
as Lg oo. For Lg = oo, the chiral symmetry is expected to be exact even at finite lattice spacing. 
Therefore, there is no need for fine tuning. DWF offer the oportunity to separate the continuum limit, 
o ^ 0, from the chiral limit, Lg —>■ oo. 

First Monte Carlo simulations 01 for AT = 1 SU{2) SYM with DWF using the HMDR l39l . 
indicate the formation of a gluino condensate which is sustained at the chiral limit. The condensate is 
non-zero even for small volume and small supersymmetry breaking mass [ 1 1 1. DWF have two welcomed 



properties that Wilson fermions do not have: There is no need for fine tuning. The second one has to do 
with the sign of the Pfaffian which is positive while in the case of Wilson fermions this is not assured ^. 

For the other hand, DWF introduce two extra parameters: Lg and uiq. These two parameters 
together with the four dimensional mass, m/, control the effective fermion mass ?7ie//- In the free 
theory one find, 

me// = mo(2 - ?Tio)[m/ + (1 - mo)'^"] , (8) 

with < niQ < 2. The value of mo = 1 is optimal because finite Lg effects do not contribute to 
rueff- In the interactive theory, one would not expect such optimal value, due to the fact that mo will 
fluctuate. Then the goal would be to have Lg large enough to have the second term of Eq. (ISj small, in 
order to simulate at reasonably small masses and extrapolate to the chiral limit, m/ — > and Lg —>■ oo 
|35|. There is a sort of "fine tuning" in order to achive the better value in Eq. DWF are much 
more expensive than Wilson fermions from the computational point of view. Wilson fermions are less 
expensive for SYM than for QCD. 



3.3 The supersymmetric Ward-Takahashi identity (WTi) 

The numerical simulations of the WTi 1 8 1 have been performed in order to determine non-perturbatively 
a substracted gluino mass and the mixing coefficients of the supercurrent. Thus, the supersymmetric 
WTi reads 0, 

{0{y)V^S^{x)) + ZTZs'{0{y)V^T^{x)) = m^Z^ , (9) 

which can be computed at fixed (3 and k. By choosing two elements of the 4x4 matrices, having 
previously choosing the operator insertion 0{y) in Eq. Q, a system of two equations can be solved for 
ZtZ^^ and muZg^ |8|. The results show a restoration of supersymmetry in the continuum limit up to 
0{a) effects. The vanishing gluino mass, extrapolated when determine rriRZg^, occurs at a value of the 
hopping parameter in agreement to the one calculated using the chiral symmetry breaking f6l. 

The supersymmetric WTi has been also studied in lattice perturbation theory lllOii40J . The idea 
is to use an independent way to determine the value of Zt/Zs, which is equal to Zt/Zs = Zxi-ioop' 
in lattice perturbation theory ^. The starting point of the calculation is the renormalized supersymmetric 
WTi on the lattice, introduced in Ref. flOl . 

Zs(OV,.S„(rc)> + Zt{OV,,T^{x)) - 2(m„ - m)Z-'{Ox''(i)> + 

^-(^l^-> - ^-("ill^-* - ^-<°|^l^-> ^^'-AOB,} =0, (,0) 

where ^(x) is a localized transformation parameter, = — ^Tr {Gpr{x) a pT-'jfj,X{x)), is a local 

definition of the lattice supercurrent which mixes with T^{x) = —'^Tr{Gpu{x)^v\{x)). x{x) = 
-^Tr {Gpr{x) a prX{x)) is the gluino mass term, is the symmetric lattice derivative, Gpr{x) is the 
clover plaquette operator and ap^ = |[7p,7r]- Then, {-^\(.=o), {O ^^k=o) and {O ^^k=o) are 
the contact terms, gauge fixing terms and Faddeev-Popov terms, respectively (we do not report them here, 
see Ref. 1 10 1). Notice that in Eq. (flOt these terms are also renormalized. This is due to the fact that their 
one-loop corrections are not just multiples of the corresponding tree-level values. J2j ^Bj (OBj) repre- 
sent the operator mixing, not only with non-gauge invariant operators (in the case the operator insertion 
O is non-gauge invariant), but also with gauge invariant operators which do not vanish in the off-shell 
regime but vanish in the on-shell one. In principle, one would require a complete list of them, or, as 
in our case, a sub-list of operators whose contributions are different from zero to the renormalization 
constant we are interested on (Zt in this case). 

'in Wilson fermions tliere are no serious sign problems when k < kc f7||8| 
^That is because, Zs = 1 + 0{gl) and Zt = 0{gl). 



In the supersymmetric limit, the renormalized gluino mass is zero, so the third term of the first 
line in Eq. (flOl vanish and we leave with a simple expression. From now on, when we refer to Eq. (flOt 
we will assume this term equal to zero. Considering now each matrix element in Eq. (flOl with O a non- 
gauge invariant operator given by O := A'l,{y) A"(z). Each matrix element in Eq. dlOt is proportional to 
each element of the Gamma-matrix base T = {1, 75, 7a, 757a, Uap}- To determine Zt one needs the 
projections over 7^ and 7a 75- 

In Fourier transformation (FT), we choose p as the outcoming momentum for the gluon field 
and q the incoming momentum for the fermion field A. To calculate Zt one should pick up from each 
matrix element of Eq. (fTOb those terms which contains the same Lorentz structure as and T^, to tree- 
level. Those operators which do not contain the same tree-level Lorentz structure as and do not 
enter in the determination of Zt- 

The renormalization constants as well as the operators, can be written as a power of go flOll . 

7 _ 7(0) , 27(2) 

^operator — ^operator "r iJO '^operator "r ' ' ' > 

{Operator) = {Operator)^^^ + gQ{Operator)^'^^ + • • • , (H) 

where {Operator)^'^\ is the 1-loop correction while {Operator)^^\ is the tree-level value. To get Zt 
we need to distinguish the tree-level values of and and for that reason we require general external 
momenta, p and q. 

Substituting Eq. (fTTl into Eq. (fTOl . up to order gg, and using the projections over 70, and 7a 75 we 
obtain, respectively fTOl . 

-tr(7a(V^S'^)^^^^p) + zP2i{paPu - Paqu - P^Sau + P • qSau) + 

Z^h{paPu - Paqu - p'^Sau + P ■ qSau) + ^tl"(7a (^||^ k=o)imp) + 

Z^Jl2i{paqu - P ■ qSau +p'^Sau) - Z^^l2ipaPu - \^(^a{^^\i;=o)arnp) 

+ ^zg^)tr(7aOi?,y°) = (12) 



and 



^tr(7a75(V^S'/,)^^p) + Zfhippq^e^apa - Z^l2ippq„euapa + 

f<^»^5i-^k=0)arnp) " 4^(7^75 k=o):,4) + 

^Zg)tr(7a750i?i>^°^=0. (13) 



More explicitly, the non-trivial part is the computation of the one-loop correction of the projections over 
Tr = {7a, 7a75} (Off-Shell regime), for tr(r,(V^5^)i'4) diagrams), tr(r,(|i^|5=o)?4) (12 dia- 
grams) and tr(rr(^|^|^=o)i^p) (4 diagrams), for a total of 28 Feynman diagrams for the determination 
of z!j?^ (see Ref. QOl for details). 

The determination of z!^^ require in principle the knowledge of which operators in Z^hT{rrOBj)^^^ 

should be included. One can demonstrate IIOI that the calculation of z!^^ does not involve any of this 
operators. Another point here is that, even in the continuum limit, Lorentz breaking terms appears in 



Eqs. (I12ll3t . but imposing the on-shell condition on the gluino mass, (once the value of Z^^ has been 
determined) they dissapear and the continuum W 
which gives good agreement with the one in 0. 



(2) 

determined) they dissapear and the continuum WTi is recovered. The final value is, Zrp \ i-ioop = 0.332, 



4. OUTLOOK 



A big effort has been made in order to describe supersymmetry on the lattice. Wilson fermions have been 
used in realistic computations with nice results. Improved chiral fermions results are starting too. New 
interesting proposals, l[T2l . are waiting for numerical applications. 
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